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Abstract: We present estimates of possible data on spin-dependent asymmetries in 
inclusive scattering of high energy polarized electrons by high energy polarized pro- 
tons at HERA with their statistical errors and discuss systematic errors. We show 
that these data will provide important information on the low-x behavior of the po- 
larized structure function gx, and will reduce the uncertainty in the determination 
of the first moment of the polarized gluon distribution Ag(x, Q 2 ) obtained from 
the QCD analysis of g\ in NLO. Furthermore, using asymmetries for di-jet events 
from a polarized HERA would substantially reduce the uncertainty in the shape of 
Ag(x, Q 2 ). Use of the information on Ag(x, Q 2 ) from the di-jet analysis in conjunc- 
tion with the NLO QCD analysis of g\ will provide an accurate determination of 
Ag(x, Q 2 ) and its first moment. 



1 Introduction 



Measurements of nucleon structure functions by lepton-nucleon inclusive electromagnetic deep 
inelastic scattering (DIS) were of fundamental importance in studying nucleon structure and 
have provided crucial information for the development of perturbative QCD (pQCD). The his- 
tory of these experiments over the past 40 years has shown that important new information has 
been obtained when measurements were extended to new kinematic regions. In the mid-1950's 
at Stanford, measurements of elastic electron-proton scattering were extended to a higher Q 2 
range of 1 (GeV/c) 2 and for the first time it was observed that the proton has a finite size [[IJ . In 
the late 1960's at SLAC the extension of measurements of inelastic inclusive electron scattering 
to the deep inelastic region of Q 2 > 1 (GeV/c) 2 led to the discovery of the parton substructure 
of the proton [Q] . This substructure was further studied in muon-nucleon deep inelastic scat- 
tering ||. The extension of the data to lower x and higher Q 2 made possible by the HERA 
ep-collider showed a surprisingly steep rise towards low x of the unpolarized structure function 
F 2 whose explanation has led to an improved understanding of perturbative QCD. Further- 
more, these data as well as data from di-jet production |J enabled a precise determination of 
the gluon density in the nucleon. 

Polarized DIS, which measures the spin dependent structure function of the nucleon has a 
similar history. The first measurements by the Yale-SLAC collaboration || extended down to 
x = 0.1 and were consistent with the naive quark-parton model view that the nucleon spin is 
carried by its constituent quarks. These experiments were followed up at CERN by the EM 
Collaboration extending the measurements to lower x of 0.01. The results showed that the 
Ellis- Jaffe sum rule is violated, which within the quark parton model implied that the contri- 
bution of the quark spins AS to the proton spin is small. This surprising result stimulated a 
large amount of experimental and theoretical work on polarized structure functions ||. The 
next generation of experiments by the SM Collaboration at CERN |9|, [10], by the SLAC col- 



laborations [12, 11, 13], and by the HERMES collaboration at DESY [14] reduced significantly 
the statistical and systematic uncertainties on the measurements and the kinematic range was 
extended to x — 0.003. In addition, spin structure functions for the deuteron and neutron 
as well as the proton were measured, thereby providing the first experimental verification of 
the fundamental model-independent Bjorken sum-rule. On the theory side, in the QCD quark 
parton model, the breaking of the Ellis- Jaffe sum rule can be interpreted as due to a positively 
polarized gluon distribution and/or a negatively polarized strange quark contribution. Hence 
precise information on the polarized gluon distribution is important for the understanding of 
the spin structure of hadrons, and will be studied in this paper. 

Despite the recent experimental progress, the polarized structure functions are still measured 
in a limited x and Q 2 range and a large uncertainty on the first moment of g± comes from the 
unknown low- a; behavior of g\. The kinematic range accessible by present-day fixed target 
experiments is limited and therefore to obtain a significant extension of the kinematic range 
the scattering of polarized electrons by polarized protons must be studied in a collider mode such 
as HERA. While the polarized quark distributions Aq are accessible in fixed target experiments 
from measurements of gi, the polarized gluon distribution Ag is deduced from scaling violations 
of g\. However, the good quality data on g 1 in the present kinematic range allows us to 
perform complete next-to-leading order (NLO) QCD analyses of the Q 2 dependence of g x . 
Recently several such analyses have been performed |15|, [16], [L7|, [Tj| and have provided the first 
information on the first moment of the polarized gluon distribution, Jq Ag(x)dx. These analyses 
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Figure 1: Measurable x — Q 2 region with polarized HERA shown with the presently explored 
regions by CERN(SMC), SLAC and DESY(HERMES) experiments, and the kinematic limit of 
measurability at HERA. 



suggest a fairly large contribution of the gluon to the spin of the proton [[[(], |1£| . Nevertheless 
the uncertainties are very large and the shape of Ag(x) is poorly constrained. This analysis 
would benefit greatly from additional data over a wider kinematic range. 

The HERA collider at DESY provides 27.5 GeV electrons and 820 GeV protons. At present 
the proton beam at HERA is unpolarized whereas the electron/positron beam has a natural 
transverse polarization, which can be rotated into a longitudinal polarization in the interaction 
regions. Once the protons in the HERA ring are polarized it will be possible to measure 
polarization asymmetries and thus explore spin dependent structure functions in the full HERA 
kinematic range. The possible extension in the measured x—Q 2 range for spin structure function 
measurements is shown in Fig. [I]. The scientific and technical problems involved in obtaining 
polarized protons in the HERA ring have been discussed in p0| . The measurability of the 
polarized structure functions at HERA has first been discussed in ETI, 1221, 123] . 



In this paper we study the potential impact of a polarized HERA on the knowledge of the 
polarized structure function and the polarized gluon distribution in the proton. We first show 
that g\ data from polarized colliding beam experiments at HERA would provide information 
on the low-x behavior of gi and through QCD analyses of gi would substantially improve the 
determination of the first moment of the polarized gluon distribution. We present for the first 
time an analysis which includes Ag(x, Q 2 ) measurements from di-jet asymmetries together with 
the QCD analysis of g\ data from fixed target experiments and from HERA. The determination 
of the polarized gluon distribution Ag(x, Q 2 ) and its first moment improves substantially when 
using such a combined analysis. 
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2 Current status of polarized parton distributions 



The structure function g± is related to the polarized quark and gluon distributions through 

- 1 dy r 



gi {x,t) = \{e 2 ) 



x y 



a a (t))AE(y,t) 



a s (t))Ag(y, t) + Cf (J, «.(*)) Ag NS (y, t) 



(1) 



where (e 2 ) = n^ 1 1]fc=i e fc, £ — m (Q 2 /A 2 ), AE and Ag NS are the singlet and non-singlet polar- 
ized quark distributions: 



AEOM) =^2^Qi{x,t), Aq NS (x,t) = £(e 2 /(e 2 > - l)A ft (x,t), 



i=l 



i=l 



and Cj NS 



a s(Q 2 )) and C g (a s ((5 2 )) are the quark and gluon coefficient functions. 

The x and Q 2 evolution of the polarized quark and gluon distributions is given by the 
DGLAP evolution equations [24]]. The full set of coefficient functions |2f| and splitting func- 
tions [26] has been computed up to next-to-leading order in a s . At NLO splitting functions, 
coefficient functions and parton distributions depend on the renormalization and factoriza- 
tion scheme. The scheme choice is arbitrary and parton distributions in different factorization 
schemes are related to each other by well-defined linear transformations. 

The polarized parton distributions can be obtained by fitting the experimental spin structure 
function data. The distributions are parametrized at an initial scale and are evolved using 
the DGLAP equations to values of x and Q 2 of the data, varying the parameters to fit the 
experimental values of g\ []15 , 16|, [17]. Here we follow the procedure used in ref. [ 15| . The initial 
parametrization at Q 2 = 1 GeV 2 is assumed to have the form 



Af{x,Q 2 ) = N{a f ,(3 f ,a f ) i\ s x a f{l - x) p * (1 + a f x), 



(2) 



where N(a, (3, a) is fixed for each distribution by the normalization condition, N(a, (3, a) Jq — 
x) p (l + ax)dx = 1, and Af denotes AE, Ag NS , or Ag. With this normalization the parameters 
VgyVNS, and 7/s are the first moments of the polarized gluon distribution, the non-singlet quark 
distribution (in the proton or in the neutron), and the singlet quark distribution, respectively, 
at the initial scale. In this parametrization af and /3f determine the low-x and high-x behavior 
of the parton distribution, respectively. The polynomial describes the medium- a; region and 
was included only for the parametrization of the singlet quark distribution. Some variations of 
the initial parton distributions were also tried in order to check the dependence of the results 
on the choice of the parametrization. The evolution was performed within the AB factoriza- 
tion scheme (in which C l g = — |^). The strong coupling constant a s (M|) was taken to be 
0.118 ± 0.003. The factorization scale M 2 and the renormalization scale /i 2 were taken to be 
M 2 = ki ■ Q 2 and /i 2 = k 2 ■ Q 2 , with k\ = &2 = 1. Further theoretical details of the fit pro- 
cedure and analysis can be found in refs. p"5L Following this scheme, the SMC performed 



a QCD analysis in NLO to existing polarized structure function data for proton, neutron and 
deuteron measured at SLAC, CERN and DESY 0, |, 0, 0, 0, 0, 0. We repeated this fit 
and show the results for the best fit coefficients in Table 1. The world data on g\ are compared 
to the best fit in Fig. 0. The first moment of the gluon distribution resulting from this fit 
is: r]g = 0.9 ± 0.3(exp) ± l.O(theory) at Q 2 = 1 GeV 2 . Similar results were obtained by an 
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Table 1: Results from the NLO QCD fit to all available published data compared to the results 
from a fit including also projected data from HERA. The parameters 77^5 are fixed by the octet 
hyperon (3 decay constants: g A = F + D = 1.2601 ± 0.0025 || and F/D = 0.575 ± 0.016 ||, 
(as = SF — D). The parameter r] g is the first moment of the gluon distribution. All values of 
the parameters are given at Q 2 = 1 GeV 2 . Errors reflect the statistical and systematic errors 
on the data. 




Figure 2: The NLO fit to proton g\ data. The solid lines are the fit to the data at measured 
Q 2 values, and the dashed and the dotted lines are the fit evolved to Q 2 = 1 and 10 GeV 2 
respectively. 
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analysis performed by Altarelli et al. [19|. The experimental error depends on the statisti- 
cal and systematic uncertainty in the measured data points, while the theoretical uncertainty 
comes mainly from the choice of the factorization and renormalization scales (Srj g = ±0.6), the 
initial parton distribution (5f] g = ±0.5) and the uncertainty in the value of a s (5r] g = ±0.3). 
The polarized gluon distribution is deduced mostly from the observed scaling violations in the 
intermediate and low-x regions. Clearly the uncertainty on the size of the polarized gluon dis- 
tribution is large. The existing data poorly constrain the low-x behavior of the various parton 
distributions and hence the values of a^,a g , and a^s- Moreover different functional forms for 
the initial parton distributions result in different predictions for the low-x behavior of g\(x) as 
well as different values of the first moment of the polarized parton distribution [Tj|, although 
the quality of the fit remains practically the same. 

The behavior of g\ at low x is still unknown. Unpolarized structure function data measured 
at HERA showed that F<i rises at low x which is in agreement with pQCD |27j, and for Q 2 > 
1 GeV 2 does not follow a Regge behavior. For polarized DIS the low-x information is crucial 



to obtain a good determination of the moments of g\ [II], |TSl [19|]. In the past the extrapolation 
of gi from the measured region down to x — 0, required to obtain its first moment, was done 



assuming a Regge behavior of the structure function, which implies |28| g\ oc x a as x — ► 



with < a < 0.5, i.e. a valence-like behavior of g\. This behavior disagrees with that predicted 
by QCD. A comparison between the Regge and QCD extrapolations from a recent analysis by 
the SMC |L0| is shown in Fig. |[ The low-x behavior has significant consequences on the first 
moment of g±. While the integral in the measured x-range of 0.003 < x < 0.7 amounts to 
0.139 ± 0.006(stat.) ± 0.008(syst.) ± 0.006(evol.) at Q 2 = 10 GeV 2 (the last error is due to the 
uncertainty in the Q 2 evolution), the contribution from the unmeasured low-x range amounts to 
0.002 or —0.011 depending on whether a Regge or a pQCD prescription is used to extrapolate 
into the low-x region. This provides the largest uncertainty on the first moment of g±. Hence, 
an experimental measurement of the low-x behavior of gi would lead to significant insight on 
the validity of QCD. 



3 Measurement of gi(x,Q 2 ) at HERA and its impact 

Presently all measurements of polarized structure functions are made by deep inelastic lepton- 
nucleon scattering on fixed targets. To evaluate the impact of measurements in the extended 
kinematic range accessible at HERA, a detailed study has been performed assuming an inte- 
grated luminosity of £ = 500 pb _1 . After the planned luminosity upgrade of HERA |31[] this 
would be achievable in about three years of data-taking. Beam polarizations were assumed to 
be 0.7 for electrons and protons. Kinematic cuts were applied based on standard analyses of 
data from the present detectors at the HERA collider |32[ [33|] : 



the angleQof the scattered electron, 9' e , is required to be < 177°; 
the inelasticity y is limited to 0.01 < y < 0.9; 

the energy of the scattered electron, E' e , is required to be > 5 GeV; 
data were taken with Q 2 > 1 GeV 2 . 



1 Angles are denned with respect to the direction of the proton beam. 
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Figure 3: The statistical uncertainty on the structure function g\ measurable at HERA for 
Q 2 = 10 GeV 2 with an integrated luminosity £ = 500 pb _1 and a comparison of various low-x 
behaviors of gf. Also shown are the SMC measurements. Starting from the measured values 
of g\ by SMC, the possible values for gf at low-x are shown for three scenarios: valence-like 
Regge behavior g\{x) oc x~ Q with a = 0, a strong powerlike positive rise g\{x) oc l/(x • (lnx) 2 ), 
and the pQCD prediction based on the present best fit. A detailed comparison in the low-x 
region of the SMC data is shown in the inset. 



A new parametrization [|Kj for the unpolarized structure function F2, which includes recent 
data from NMC, E665, HI and ZEUS was used. The SLAC parametrization |34| was used for 
the ratio R of longitudinal and transverse 7*p cross sections. Results for the expected number 
of events for different x — Q 2 bins are listed in Table 2, together with the average y, the depo- 
larization factor D = (y{2 — y))/{y 2 + 2(1 — y)(l + R)), and the expected statistical uncertainty 
on the measured asymmetry 5A m . Standard deep inelastic scattering formulae ^] were 
used to calculate the kinematics and statistical uncertainties on the measured asymmetries. 
Using the best-fit parton distributions obtained from the QCD analysis of the data discussed in 
Sect. 2, predictions were made for gi(x,Q 2 ) values and hence for cross section asymmetries in 
the HERA range. Systematic uncertainties were not included and will be discussed at the end 
of this section. As shown in Fig. [1], measurements of g\ at HERA will extend the x region down 
to x = 5.6 x 10~ 5 . In Fig. |3| we show the expected accuracy of the determination of g\ at HERA 
with the maximal extent of variation in the low-x behavior compatible with the requirement of 
integrability of g\ (which implies that at small x g\ can rise at most as l/(x ln a x) with a > 1). 
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Table 2: The kinematic variables x,Q 2 ,y, and D, with the number of events expected iVtotab 
and the statistical uncertainty in the measured asymmetry 5A m assuming an integrated lumi- 
nosity C = 500 pb _1 and proton and electron beam polarizations P p = P c = 0.7. 

For the HERA measurements only one data point is shown at each x value, chosen to be the 
point with the smallest statistical errors which corresponds to the largest depolarization factor 
(see Table 2). The figure also shows that a polarized HERA can clearly distinguish between 
the different prescriptions for the g\ behaviors in the low-x region. 

The measurable points with their estimated errors for HERA experiments are shown in 
Fig. |] for eight accessible x-bins (out of the nine bins listed in Table 2). They are predicted on 
the basis of the NLO fit of Table 1, represented by the solid line. We also consider two limiting 
scenarios: a) the first moment of the gluon distributions is fixed to be at Q 2 = 1 GeV 2 
(minimal gluon: dashed lines in Fig. |]) and b) the first moment of the singlet quark density 
was fixed to rj q = a§ at the same reference scale (maximal gluon: dotted lines in Fig. Of 
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Figure 4: Predicted values for g\ from NLO fits together with estimated statistical uncertainties 
for a future polarized DIS experiment with an integrated luminosity £=500 pb -1 and with 
beam polarizations P e = P p = 0.7. The solid lines are predictions based on fits to present data 
extended into the HERA kinematic range, while the dotted and dashed lines are the predictions 
with a maximal and minimal gluon, respectively (see text). 



course still wider deviations from these fits are possible since the low-x behavior of the current 
best fit is only very loosely constrained due to the lack of direct experimental information at 
low x. Furthermore, higher order corrections beyond NLO may turn out to be important at 
very low x |36|, |37| . 

The polarized gluon distribution is deduced from the scaling violations of g±. Therefore a 
reasonably wide kinematic coverage is required in order to achieve such a determination with 
satisfactory accuracy. In particular, most of the present-day parametrizations of polarized 
parton distributions ]T5|, [IB], [T7| indicate that Ag(x, Q 2 ) is peaked in the unmeasured \ow-x 



region. Therefore data in this region with a large span in Q 2 , such as those obtainable at 
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HERA, would substantially improve the determination of Ag(x,Q 2 ). 

To assess the impact of these data we have repeated the QCD analysis described in Sect. 
and included the projected HERA data. The values of g\ at HERA predicted by the QCD 
analysis were randomized within their estimated statistical errors by assuming a Gaussian 
distribution. The values of the best-fit parameters (Table 1) of course remain the same within 
the errors, but the calculated statistical errors on the projected measurements, used in the 
QCD analysis provide an estimate on the reduction of the experimental uncertainty in the 
fitted parameters. In Table [1] we compare the result of the fit with the projected HERA data 
to the analysis of present-day data only. Note in particular the sizable improvement in the 
determination of the first moment of Ag, from the present error 5(rj g ) = ±0.3 to ±0.2 with 
HERA data. A similar improvement is observed for the uncertainty in a g , which is determined 
by the \ow-x data. The theoretical uncertainty on the first moment of Ag as of today is large. 
The major contributing sources are factorization and renormalization scales, the initial parton 
distributions and the uncertainty in the value of a s . At least the first two are direct consequences 
of the fact that g\ structure function measurements are available from a rather small x and 
Q 2 -range, namely 0.003 < x < 0.8 and 1 < Q 2 < 100 GeV 2 . With HERA the measured x and 
<5 2 -range would increase substantially, 10~ 5 < x < 0.8 and 1 < Q 2 < 10 4 GeV 2 . Consequently 
the principal sources of theoretical uncertainty are expected to reduce significantly. Out of all 
the systematic uncertainties we tried to study the most dominant one, namely the factorization 
and renormalization scale uncertainty. This uncertainty is related to the lack of availability of 
higher than next-to-leading order corrections and is usually estimated by changing the scale 
factors, k\, k% (see Section 2), in both directions (up and down), e.g. by a factor of two Jl9[ or 
four [[U| . The difference between the central value if g st and the one obtained from the fit which 
gives the largest deviation from rf g est is assigned as the uncertainty on r] g due to the uncertainty 
in scales, 5 sca i e (r] g ). 

If we follow this procedure using the present day data, the S sca _i e (r] g ) « ±0.6. In order to 
estimate the improvement one can expect in this uncertainty from future HERA data, this 
procedure was repeated with the projected data in Table 2. The QCD analysis fails to converge 
with the scale factors changed by factors of 4 or even 2, indicating that with HERA data, 
this uncertainty is constrained already. We repeated this procedure further by changing the 
scale factors by smaller amounts ([k\, k^\ = 0.5, 0.8, 1.2, 1.5) with the requirement that the fit 
procedure should converge and result in fits with comparable x 2 -confidence level. The largest 
difference in value of g g obtained from such a procedure compared to the one in our present 
best fit was £™k A = ±0.15. This suggests a possible improvement in this source of uncertainty 
by a factor of ~ 4. This improvement is indicative of the reduction one could achieve in the 
total theoretical uncertainty in t] g from HERA data. 

Systematic errors are discussed in [^, |23|]. The principal types of systematic errors associ- 
ated with spin dependent asymmetry measurements are normalization errors and false asymme- 
tries. The normalization errors are mainly due to uncertainties in the electron and proton beam 
polarizations P e and P p . These lead to a change in the magnitude of the measured asymmetry 
by an amount which is small compared to the statistical error, and hence they are important 
primarily when evaluating the first moment of g\ (x). 

Measurements of the electron polarization P e at HERA by Compton scattering from polar- 
ized laser light have achieved [|14], |38| a relative uncertainty of 5%, and plans are being made 
to improve it to the 2% level which has already been achieved at SLAC |j39| . Absolute mea- 
surement of the polarization of the high energy proton beam presents a new challenge and is 



9 



presently under investigation at RHIC and at HERA. It is reasonable to expect that an ac- 
curacy of 5% can be achieved [[40], (4l|, |42|, f43[. The polarization errors of 2% (5%) translate 



directly into a relative error of 2% (5%) on the structure function gi. 

The other type of systematic error, which is due to false asymmetries, arises from time 
variations in detector efficiencies between the measurements for the parallel and antiparallel 
spin conditions, and the correlation of proton polarization with the beam intensity or the bunch 
crossing angle. With regard to the time variation of detector efficiencies, the interaction region 
for proton and electron bunches should be the same for different proton spin orientations, thus 
keeping the detector acceptance constant. In the polarized HERA collider we expect to run 
in a mode where a sequence of about 10 proton bunches will have a given polarization, and 
the next sequence will have the opposite polarization. Since the time interval between the 
successive bunch crossings is 96 ns, the asymmetry measurements will be done within 2 fis and 
for such a time scale the detector efficiencies should be constant. A possible correlation between 
proton bunch intensity and polarization can lead to a false asymmetry, which can be avoided 
by repeating the measurements after rotation of all the proton spins. Small contributions 
to false asymmetries may arise due to the correlations with the proton spin direction of the 
micro-structure and polarization distribution within a bunch which can also be eliminated by 
the proton spin rotation. A further possible source of false asymmetries comes from potential 
asymmetries of background contamination in the ep sample, namely interactions of the beams 
with beam-gas. Present analyses at HERA have achieved data samples with a purity of better 
than 99%. Furthermore, residual asymmetries caused by such background interactions can be 
efficiently monitored by using so-called pilot bunches, i.e. bunches which pass the interaction 
region and have no collision partner in the other beam. Therefore this source of systematics 
can be kept smaller than 10~ 4 . Hence we believe that false asymmetries at polarized HERA 
can be controlled at the required level of 10~ 4 or less. Similar requirements are necessary for 
the RHIC-SPIN program and are being actively studied. More quantitative estimates can not 
be made at this time, since they depend on the polarized collider parameters and operating 
conditions which are being studied at DESY. 

We conclude that introduction of systematic errors will not change significantly our con- 
clusions on the impact of future HERA g\ data. This conclusion also applies to radiative 
corrections, which have been studied for the HERA kinematics in [El] and we checked, using 



the program HECTOR 1.11 [45], that they are well under control. 



4 A direct measurement of Ag(x, Q 2 ) from di-jet events 
at HERA 

From the discussion in Sect. 3 it appears that the polarized gluon distribution Ag(x, Q 2 ) is 
not as well constrained by the Q 2 dependence of the structure function as is the unpolarized 
gluon density. In fact, while the first moment of the polarized gluon density can be deduced 
from scaling violations (including future HERA data) within 20%, the shape of Ag(x) is not 
constrained. Therefore measurements which are directly sensitive to Ag(x,Q 2 ) are of great 
interest and are expected to add important information to the determination of not only the first 
moment of Ag but also of its shape. It has been shown by HI at HERA || that in leading order 
(LO) the unpolarized gluon momentum distribution, g(x), can be accessed in a direct way by 
measuring the cross section of di-jet production in DIS in a region of Q 2 = 5-50 GeV 2 with both 
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jets measured in the central detector (the proton remnant jet is not counted), and for jets with a 
transverse energy E T larger than 5 GeV. In leading order two processes contribute to the di-jet 
events, Photon-Gluon Fusion (PGF) and QCD-Compton scattering (QCDC), both leading to 
two jets back to back in azimuth in the hadronic center-of-mass system. The Feynman diagrams 




Photon Gluon Fusion QCD — Connpton 



Figure 5: Feynman diagrams for the processes contributing to the di-jet cross section at LO: 
the Photon-Gluon Fusion (PGF) process (left) and the QCD-Compton (QCDC) process (right). 

for the two processes are shown in Fig. |5|. The PGF cross section is proportional to a s g(x). The 
QCDC cross section, however, is proportional to a s q(x) and constitutes the background. The 
kinematic parameters to describe the PGF process are the momentum fraction of the proton 
carried by the gluon x g , the four-momentum transfer Q 2 and the invariant mass of the two jets 
s. They are related to the Bjorken-x, by: 

For s > 100 GeV 2 , x g is larger than x by about an order of magnitude and the accessible range 
at HERA is about 0.002 < x g < 0.2. 

The feasibility of using this kind of measurement in the case of polarized lepton-nucleon 
scattering to extract Ag(x,Q 2 ) from di-jet asymmetries has been discussed in [47|, f48|. A new 



study of such a measurement at HERA has been performed and is detailed in ||49|| . Recently 
NLO corrections were calculated for the polarized jet rates and found to be small |5{J. In 
this paper we summarize the main results and use the projected data to assess the sensitivity 
to the shape of Ag(x,Q 2 ). In LO the total unpolarized di-jet cross section is the sum of the 
contributions from PGF processes and QCDC scattering, and can be written as: 

^di-jet = + = Ag + Bq (3) 

where g and q are the gluon and quark densities. A and B can be calculated in perturbative 
QCD. Similarly, for the polarized case we can write: 

A(7 di _ jct = a\i }ct - o-ff_jet = a Ag + b Aq (4) 

with Ag and Aq being the polarized gluon and quark densities. The left arrow in the subscript 
denotes the polarization of the incoming lepton with respect to the direction of its momentum. 
The right arrow stands for the polarization of the proton being anti-parallel or parallel to the 
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polarization of the incoming lepton. The di-jet asymmetry is sensitive to Ag/g, especially at 
low x where the PGF cross section dominates. 

^. jet = ^i = ^^^3L +B ^i(l-^), (5) 
^C"di-jct g 2cr di _j ct q z o" di _j et 

with A = a/A and B = b/B. The experimentally accessible asymmetry A meas is related to 
Adi-j e t by 

^meas = _|_ jyff = ^cPpD A di— jet (6) 

where D is the depolarization factor and A meas is the measured asymmetry. The quantities N*^ 
(N^) are the total number of observed di-jet events (N^ = nIq F + Nq CBC ) for antiparallel 
(parallel) proton and electron spins. 

A Monte Carlo study was performed for a polarized HERA on the expected size and uncer- 
tainty of the measured di-jet asymmetries and the contributions from signal and background 
processes. PEPSI 6 . 5 J51J was used, which is a full LO lepton-nucleon scattering Monte Carlo 
program for unpolarized and polarized interactions, including fragmentation, and unpolarized 
parton showers to simulate higher order effects. On the parton level the results were cross 



checked using the Monte Carlo cross section integration program [p2[ . Detector smearing and 
acceptance effects were included based on the performance of the present HI detector. As input 
for the unpolarized parton distributions we used the LO GRV parametrizations |k| and for the 
polarized distributions various published parametrizations, as well as the result for Ag of the 
NLO fit described in the previous section were used. All these distributions are in agreement 
with the present data. Monte Carlo events were generated for an integrated luminosity of 500 
pb _1 and beam polarizations of 0.7 for electron and proton beams. Jets were defined by the 
cone jet algorithm with a cone radius of 1 in the pseudorapidity rj and azimuthal angle <fi plane. 
Di-jet events were selected by the following criteria: 

• two jets with E T > 5 GeV; 

• 5 GeV 2 < Q 2 < 100 GeV 2 ; 

• 0.3 < y < 0.85; 

• s > 400 (100) GeV 2 . 

The invariant mass of the two jets is reconstructed from their four- momenta Pi(2)' s = (Pi + 
P2) 2 - The cut s > 400 GeV 2 was chosen to keep higher order corrections small ||50|| . However, 



the data with 100 < s < 400 GeV 2 were also analyzed. The asymmetries were found to be as 
large el. S 9 few percent and are listed in Table Ha. 

The two lowest x 5 -bins in the table are obtained from the lower s data. In this table the 

asymmetry corrected for QCDC background (A COTr = A meas Q< ^u + jv?T CDC ) i s a ^ so presented. 

Note that there is a significant QCDC contribution only for the two highest x 5 -bins. Fig. 0] 
shows the projected values of Ag{x) / g{x) and xAg(x) as a function of x g for £ = 500 pb~ 
for different Ag(x, Q 2 ) distributions: the polarized gluon density obtained from the NLO fit to 
gi data by the SMC, the gluon set C by Gehrmann and Stirling JHJ and an instanton-induced 



gluon proposed by Kochelev |5^| . The results are presented at Q 2 = 20 GeV 2 which is very 
close to the average Q 2 of the simulated data. 
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0.006 




0.034 
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-0.029 


0.007 


0.084 
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0.012 




0.084 


-0.041 


0.013 


-0.035 


0.013 


0.207 


-0.023 


-0.047 


0.027 




0.207 


-0.091 


0.027 


-0.013 


0.031 



Table 3: a) Expected measured asymmetries from di-jet events for six x 9 -bins with a polarized 
HERA and background corrected asymmetries with statistical errors corresponding to an inte- 
grated luminosity of 500 pb _1 . b) The background corrected asymmetries for two different Q 2 
ranges. 




X X 



Figure 6: The sensitivity to Ag(x)/g(x) (a) and xAg(x) (b) from a measurement of di-jet 
events with a polarized HERA at Q 2 = 20 GeV 2 assuming an integrated luminosity of 500 pb _1 . 
The error bars show the expected statistical error. An estimate of systematic uncertainties 
is reflected in the shaded error band. As input for the parton distributions three different 
parametrizations were used (see text). 



The shaded error band shows an estimate of the systematic uncertainty considering the 
following error sources: an uncertainty of 2% on the calibration of the hadronic energy scale, 
the error on the total unpolarized di-jet cross section of 2% and on the unpolarized gluon 
density g(x) of 5% [55]. The uncertainty on the ratio of the polarized and unpolarized quark 
densities, Aq/q, was considered to be 10% and the error on the polarization measurements for 
electron and proton beams was taken to be 5%. Both statistical and systematic errors were 



calculated |49j using as input the gluon set A of Gehrmann and Stirling |17|, which resembles 
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the gluon distribution extracted from the NLO fit to g\, and were adjusted to the different 
parametrizations used here. 

It can be seen that the measurement is sensitive to the different gluon distributions, and 
that its shape can be determined with six bins in x. Using these projected measurements the 
statistical error of the first moment of Ag in the x-range of 0.0015-0.32 is: 5(/ ° 3 2 15 Ag(x)dx) = 
0.2. 

A further study was performed using additional information on the Q 2 dependence. For 
this study we extended the Q 2 range down to 2 GeV 2 and divided the data into two bins with: 
2 < Q 2 < 10 GeV 2 and 10 < Q 2 < 100 GeV 2 . The expected sensitivity on Ag(x,Q 2 ), reflected 
by the expected asymmetries and their errors, is shown in Table The mean Q 2 values for 
the two bins are 4.5 and 30 GeV 2 , respectively. The effect of the Q 2 evolution of the gluon can 
be seen, but does not appear to be very large. 

5 Combined QCD analysis of gi data and Ag from di-jet 
asymmetries 

The measurement of gi(x,Q 2 ) and di-jet asymmetries provide complementary information on 
the polarized gluon distribution. As argued in Sections 3 and 4 the first moment of the polarized 
gluon distribution can be obtained from the QCD analysis in NLO of g\ data while di-jet 
asymmetries will provide a direct measurement of the polarized gluon density Ag(x, Q 2 ) in the 
range 0.0015 < x < 0.32. Hence, it is interesting to study the impact of a combined analysis of 
gi and di-jet asymmetry data on the determination of Ag. 

In principle, such an analysis should be carried out using the measured di-jet asymmetries 
and g\ data in a self-consistent procedure in NLO. Polarized quark and gluon densities should 
be fitted simultaneously to both data sets, in order to reproduce the scaling violations in g± 
and separate the PGF and QCDC contributions to the jet asymmetry cross sections. The tools 
to perform such analysis are not available yet. However, the effect of the combined analysis on 
the uncertainty in Ag can be assessed by including the values of Ag(x, Q 2 ) obtained from the 
di-jet analysis in the NLO QCD fit of the gi data as an extra constraint. 

The values of Ag(x, Q 2 ) in the six x-bins, as obtained from the analysis of di-jet asymmetries, 
were considered as data. Statistical errors were obtained from the Monte Carlo simulation 
discussed in Sect. 4. The values used for Ag(x, Q 2 ) were taken from the best fit of present-day 
gi data described in Sect. 2. Like the projected g\ data for the HERA kinematic region, they 
were randomized around the calculated value within one standard deviation of the statistical 
error. As mentioned in Section 3, including the projected g\ data in the NLO analysis, reduced 
the error on the first moment of Ag by about 30%, yielding r] g = 0.9 ± 0.2, down from S(r] g ) = 
0.3. The QCD analysis of the present and the projected g\ data was repeated including also 
Ag(x,Q 2 ) values from the di-jet asymmetries in the fit. While the best fit parameters are 
obviously about the same as those obtained for the NLO fit of the gi data, the uncertainty on the 
first moment of the polarized gluon distribution is reduced by 50%, resulting in r\ g = 1.0 ± 0.1. 
The polarized gluon distribution from the best fit of g± data, and projected HERA data, evolved 
to 20 GeV 2 is shown in Fig. ^ along with the projected data of Ag(x, Q 2 ) used in the fit. The 
effect of the di-jet data alone on the uncertainty in r] g is similar to that of the additional 
gi data which will be measured at HERA. The combined fit of present-day g\ data with the 



14 




-5 -4 -3 -2 -1 -5 -4 -3 -2 -1 

10 10 10 10 10 1 10 10 10 10 10 1 



X X 

Figure 7: The polarized gluon distribtion at Q 2 = 20 GeV 2 obtained from a 14 parameter fit 
(Eq. 3), dashed line, including only existing g\ data (a), and including existing and projected 
gi and Ag(x) data (b). For comparison also the results of the 10 parameter fit to existing and 
projected data is shown, solid line, as well as the corresponding Ag values with their statistical 
errors, as they could be obtained from a di-jet analysis. 



Data used 


S(fAg) 


QCD Anal, of published g 1 today 
(= "A") 


±0.3 


"A" + Projected HERA g\ 


±0.2 


"A" + Projected HERA A^-jets 


±0.2 


"A" + Projected HERA g\ and A^-jets 


±0.1 



Table 4: The expected uncertainty in the determination of the first moment of the gluon 
distribution at Q 2 = 1 GeV 2 with inclusive and di-jet data from polarized HERA with C = 500 
pb- 1 . 

Ag(x,Q 2 ) from di-jet asymmetries yields r\ g = 1.0 ± 0.2. The results for the first moments of 
the polarized gluon distributions obtained from the fit of g\ data with projected g\ and Ag 
data are summarized in Table 4. 

While the first moment of the polarized gluon distribution can be deduced from scaling 
violations in gi, the QCD analysis does not constrain its shape. This can be demonstrated by 
introducing an additional term, bf^Jx into the functional form of the parton distributions, 

Af(x,Q 2 ) = N(a f ,f3 f ,a f ) Vf x a f{l - xfs{\ + a f x + b f y/x). (7) 
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Using this form in the QCD analysis of the gi data results in a fit with a node in Ag(x, Q 2 ), 
which is not allowed by the more restricted form of Eq. 2. The gluon distribution obtained 
from a fit with this functional form is shown in Fig. [7] where it is also compared with the 
gluon distribution obtained with the standard form (Eq. 2). While the first moment of Ag at 
Q 2 = 1 GeV 2 is 0.9 ± 0.3, the same value as obtained in the standard fit, the shapes of the 
gluon distributions are different. 

On the other hand, di-jet asymmetries provide good information on the shape of the po- 
larized gluon distribution in the measured region. The effect of the additional data on our 
knowledge of Ag can be demonstrated when the projected values of g-y and Ag(x,Q 2 ), (ob- 
tained from the standard, 10 parameter fit) are fitted with 14 parameters, as discussed above. 
As shown in Fig. [7|, once all the projected g\ and di-jet data are included in the fit, the po- 
larized gluon distribution obtained from the 14 parameter fit is essentially identical to that 
obtained from the 10 parameter fit in the region where Ag(x, Q 2 ) data are introduced. Obvi- 
ously, Ag(x, Q 2 ) remains unconstrained in the high x region where no Ag(x, Q 2 ) data is added 
and where the Q 2 evolution is small. One may thus conclude that the combined analysis will 
constrain the functional form of the polarized gluon distribution in addition to its first moment. 

6 Conclusions 

Inclusive polarized DIS measurements at HERA with polarized protons and with a high inte- 
grated luminosity will yield significant and unique new information on gf(x, Q 2 ) over a much 
extended range in x and Q 2 , which is of great theoretical interest. Statistical errors on the points 
will dominate. For such measurements the false asymmetries should be considerably smaller 
than the true asymmetries to be measured, and normalization errors should be controlled to 
be better than 10%. 

The present uncertainties in the determination of the first moment of the polarized gluon 
distribution are large. They will be considerably reduced by the pQCD analysis of the g\ 
structure function data from HERA. However we have shown that even further reduction in 
the uncertainty can be achieved if polarized gluon information coming from the di-jet events is 
used in the analysis in conjunction with the g\ data. 

The fact that HERA is an operating facility with a plan for substantial increase in luminosity 
and with two major operating detectors, as well as a polarized electron beam, means that only 
the high energy polarized proton beam needs to be developed. The efforts to achieve this 
goal seem well justified and data from such experiments at HERA would provide unique and 
complementary information to that from other presently proposed experiments. 
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